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We study U(l) gauge theory on a 4d non-commutative torus, where two directions are 
non-commutative. Monte Carlo simulations are performed after mapping the regularized 
theory onto a U(A'^) lattice gauge theory in d = 2. At intermediate coupling strength, 
we find a phase in which open Wilson lines acquire non-zero vacuum expectation values, 
which implies the spontaneous breakdown of translational invariance. In this phase, various 
physical quantities obey clear scaling behaviors in the continuum limit with a fixed non- 
commutativity parameter 6, which provides evidence for a possible continuum theory. In 
the weak coupling symmetric phase, the dispersion relation involves a negative IR-singular 
' term, which is responsible for the observed phase transition. 

o 

O ■ §!• Introduction 



Non-commutative (NC) geometry has been studied extensively as a modification 
of our notion of space-time at short distances. It has recently attracted much at- 
tention since gauge theories on a NC geometry have been shown to appear as a low 
energy limit of string theories with a background tensor field.'^^ NC geometry also 
appears naturally in matrix model formulations of string theory .'^J® Introducing 
non-commutativity to the space-time coordinates may change the infrared dynamics 
drastically at the quantum level due to the so-called UV/IR mixing.^' The same 
effect also poses a severe problem in the renormalization procedure within perturba- 
tion theory since a new type of IR divergences appears in non-planar diagrams. 

The finite lattice formulation^' (extending an earlier worliP) regularizes such 
divergences as well as the ordinary UV divergences. It therefore provides a non- 
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perturbative framework to establish the existence of a consistent field theory on a 
NC geometry. In ref. 6} a simple field theory — 2d U(l) gauge theory — on a 
NC plane has been studied by Monte Carlo simulation, and the existence of a finite 
continuum limit has been confirmed. The ultraviolet dynamics is described by the 
commutative 2d U((X)) gauge theory. On the other hand, Wilson loops of large area 
(with a variety of shapes) pick up a complex phase linear in the loop area, in the spirit 
of the Aharonov-Bohm effeclj^ (but their expectation values are shape-dependent, 
in contrast to the commutative plan^). 

As an interesting physical consequence of the UV/IR mixing in the case of scalar 
field theory, ref.Ej) predicted the existence of a "striped phase" for the NC Xcf)^ model 
in dimensions d> 2. In this phase, non-zero Fourier modes of the scalar field acquire 
a vacuum expectation value and break the translational invariance spontaneously. 
The existence of this new phase was fully established by Monte Carlo simulations 
Such a phase has also been observed in the 2d lattice model,1^''^ but in that case 
the continuum limit is still under investigation. An analogous phase was observed 
in simulations involving a fuzzy sphere instead of a NC plane.'^ 

In the case of 4d U(l) NC gauge theories, the one-loop calculation of the effective 
actiori^^' suggests that the perturbative vacuum is unstable against the condensa- 
tion of the Wilson line. This causes the spontaneous symmetry breaking (SSB) of 
the translational invariance, since the open Wilson lines carry non-zero momenta.!^ 
Whether a stable non-perturbative vacuum exists or not is an interesting question, 
which we have addressed in ref. [T^ from first principles using the lattice formula- 
tion.El 

§2. Lattice gauge theory on a NC geometry 

NC geometry is characterized by the commutation relation = iOf^u 

among the space-time coordinates x^, where is the non-commutativity ten- 
sor. Here we consider the 4d Euclidean space-time {fi = ,4) with the 
non-commutativity introduced only in the /U = 1,2 directions, i.e. 

©12 = -021=^, (2-1) 

and 0^1, = otherwise. Since we have two commutative directions, we may regard 
one of them as the Euclidean time. This allows us to alleviate the well-known 
problems concerning causality and unitarity. 

The lattice regularized version of gauge theory on NC geometry can be defined 
by an analog of Wilson's plaquette actiorP 

S = -p^^U^{x)i^U^{x + afi)*U^{x + ai))* i.U^{x)* + c.c. , (2-2) 

where the symbol jl represents a unit vector in the /i-direction and we have in- 
troduced the lattice spacing a. Here the star product (denoted by -k) encodes the 
non-commutativity of the space-time. 

In order to study the lattice NC theory ()2-2p by Monte Carlo simulations, it 
is crucial to reformulate it in terms of matrices.!^ In the present setup (|2-ip . with 
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two NC directions and two commutative ones, the transcription applies only to the 
NC plane whereas the commutative plane remains untouched. Let us decompose 
the four-dimensional coordinate as x = (y, z), where y = {xi,X2) and z = {x3,X4) 
represent two-dimensional coordinates in the NC and in the commutative plane, 
respectively. We use a one-to-one map between a field (p{x) on the four-dimensional 
N X N X L X L lattice and a, N x N matrix field (p{z) on a two-dimensional L x L 
lattice. This map yields the following correspondence 

ipi{y,z)i.ip2{y,z) <^ ipi{z)ip2{z) , (2-3) 
ip{y + aft, z) <^ ip{z) rl , (2-4) 

l^Y.^iy,z) ^ ^tr^(z). (2-5) 

y 

The SU(A^) matrices -T^ {fi = 1, 2), which represent a shift in a NC direction, satisfy 
the 't Hooft-Weyl algebra with the matrix size being odd. For this construc- 
tion,l3'l3 the non-commutativity parameter 9 in (j2Tp is given by 

e = -Na'^ . (2-6) 
vr 

Note that the extent in the NC directions Na goes to infinity in the continuum limit 
a ^ at fixed 6. 

Using this map, the link variables U^{x) are mapped to a x unitary matrix 
field U^{z) on the two-dimensional L x L lattice. The action (l2-2p can be rewritten 
in terms of U^{z). By performing a field redefinition V^(z) = tJ^{z)r^ for = 1,2 
and V^{z) = U^iz) otherwise, we arrive at 

S = SnC + 'S'com + 'S'mixed 5 

Snc = -NI3Zi2Y,iT {Vi{z)V2{z)Vi{z)^ V2{z)^)+ c.c. , 

z 

5com = -A^/3^tr [v3{z)Va{z + a%)V3{z + al)^ V^iz)^) + c.c. , 

z 

2 4 

5mixed = -iV/?J]J]J]tr (v^,{z)V,{z)V^{z + av)^V,{z)^) + c.c. . (2-7) 

z {1=1 u=3 

This action possesses the U(l)^ symmetry Vfj_{z) i— > e*°^V^(2:) for fj, = 1,2, which is 
related to the translational symmetry in the NC directions of the action (j2-2p . 

§3. Phase diagram 

Let us first investigate the phase structure of the lattice model (|2-7p . As an 
order parameter for the spontaneous breaking of the U(l)^ symmetry, we define a 
gauge invariant operator 
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for ^ = 1,2. This operator corresponds to the open Wilson hne carrying a momentum 
with the absolute vshi^^ 

— h — I t ^"-"^ 's.^en n , , , 

^ " ^ ' - \ -tiV for odd n . ^''■^^ 

Since the operator P^{n) with odd n carries a momentum on the cutoff scale, it 
does not couple to excitations that survive in the continuum limit. Therefore we will 
focus mainly on the even n case in what follows. 

In fig. [1] we plot (|Pi(n)|) against j3 
for n = 2. We observe that there is a 
phase, in which the order parameter be- 
comes non-zero. On the other hand, the 
quantity (|Pi(n)|) for odd n takes tiny 
values throughout the whole region of 
(5. This implies that the U(l)^ symme- 
try is broken down to (Z2)^ in this phase, 
which we refer to as the "broken phase" . 
The critical point between the broken 
phase and the weak coupling phase, de- 
noted as /3c, increases as (3c ~ A^^- 

§4. Continuum limit 

In this section we investigate 
whether it is possible to fine-tune f3 as a 
function of N in such a way that phys- 
ical quantities scale. For that purpose, 
let us next consider closed Wilson loops, 
which play an important role in commu- 
tative gauge theories as a criterion for confinement. In the present case, since we 
introduce non-commutativity only in two directions, there are three kinds of square- 
shaped Wilson loops depending on their orientations. For instance, the Wilson loops 
in the NC plane and the commutative plane are defined respectively as 

z 

Wu{n) = ^Y.^v{y3{z,n)Vi{z + na^,n)V^{z + nal,n)^Vi{z,n)^) , (4-1) 

z 

using the parallel transporter in the commutative directions 

Vu{z, n) = K(z)K(^ + az)) • • • V^,{z + (n - l)oi>) . (4-2) 

In the following we set a = 1 for = 45 as a convention, and the lattice spacing 
a for other is determined through (|2-6p with = ^ — 14.3. As a practical strategy 
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Fig. 1. The order parameter (jPi(n)|) is plot- 
ted against {3 for n = 2. The system size 
is TV = 15 (circles), N = 25 (triangles) 
and A'^ = 35 (squares). The closed (open) 
symbols represent results obtained with in- 
creasing (decreasing) /3, which show a clear 
hysteresis behavior. 
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Fig. 2. The expectation value of the Wilson 
loop in the commutative plane. 
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Fig. 3. The expectation value of the Wilson 
loop in the NC plane. 



to fine-tune /3, we first optimize the scaling behavior for the expectation value of the 
square Wilson loop W34(n) in the commutative directions. Fig. [2] shows the result. 
The horizontal axis represents the physical side length na of the loop. We observe 
a clear scaling behavior, and the scaling region extends as increases. In fact the 
optimal (3 increases with N much slower than the lower critical point (3c between the 
broken phase and the weak coupling phase, which grows as N"^. This implies that we 
remain in the broken phase in the continuum limit. In fig. [3] we plot the expectation 
value of the Wilson loop in the NC plane with the same sets of parameters as the 
ones used to obtain fig. [2j We do observe a compelling scaling behavior. We have 
also confirmed the scaling of other quantities.'^ 



In this section we discuss the dispersion relation in the symmetric phase. As 
we mentioned before, in the present setup we regard one of the commutative coordi- 
nates (say, 3:4) as the Euclidean time. From the exponential decay of the two-point 
correlation function of open Wilson line operators separated in the time direction, 
we can extract the energy of a state that couples to this operator. Similar studies 
have been done also in the case of NC scalar field theory.'^J 

Let us define the open Wilson line operator at a fixed time X4 as 



for /i = 1,2, which has a zero momentum component in the X3 direction and a non- 
zero momentum component (j3-2p in a NC direction depending on n. Then we define 
the two-point correlation function of the open Wilson lines 



§5. Dispersion relation 




Cn{r) ^ j;(P^(x4,nr •P^(x4 + T,n)) 



(5-1) 



with a separation r in the temporal direction. 



6 J. Nishimura, W. Bietenholz, Y. Susaki, J. Volkholz 



Fig. m shows the dispersion relation 
obtained from the two-point correlation 
function (|5-ip for even n. The param- 
eters {N, (3) are chosen as in the bro- 
ken phase. Namely, a is still determined 
through dMI with 61 = ^ ~ 14.3, and 
we fine-tune (3 at each N in such a way 
that the scaling behavior of the Wilson 
loops in the commutative plane is opti- 
mized. It turns out that the data points 
{E,p) for different lie to a good ap- 
proximation on a single curve 
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Fig. 4. The dispersion relation in the symmet- 
ric phase. The energy E obtained from the 
(5-2) two-point correlation function (|5-ip is plot- 

ted against the momentum p. 

with c ~ 0.1285. Due to the negative 

sign of the second term in (|5-2p . the usual Lorentz invariant (massless) dispersion 
relation is bent down. The IR singularity is regularized on the finite lattice since the 
smallest non-zero momentum (which corresponds to n = 2) is given by oc 
However, if one increases at fixed 9, the energy at the smallest non-zero momentum 
vanishes at some A^, and one enters the broken phase. Therefore we cannot take 
the continuum limit in the symmetric phase. We have also studied the dispersion 
relation in the broken phase, which reveals the existence of a Nambu-Goldstone mode 
associated with the SSB of the U(l)^ symmetry-^SJ 



§6. Summary and discussion 



To summarize, we studied four-dimensional gauge theory in NC geometry from 
first principles based on its lattice formulation. In particular we clarified the fate 
of the tachyonic instability encountered in perturbative calculations. The lattice 
formulation is suited for such a study since the IR singularity responsible for the 
instability is regularized in a gauge invariant manner, and we can trace the behavior 
of the system as the regularization is removed. This revealed the existence of a 
first order phase transition associated with the spontaneous breakdown of the U(l)^ 
symmetry, which corresponds to the translational symmetry in the NC directions. 

The dynamical extent in the NC directions — defined through the eigenvalue 
distribution of matrices — turns out to be finite in the continuum limit .'^ An 
analogous first order phase transition is found in gauge theories on fuzzy manifolds,'^ 
where the fuzzy manifolds collapse at sufficiently large couplings. The instability 
in those cases is due to the uniform condensation of a scalar field on the fuzzy 
manifold. The phenomenon that the space-time itself becomes a dynamical object 
is characteristic to gauge theories on NC geometry. This may be related to the 
dynamical compactification of extra dimensions in string theory .l^^ See also ref. llTh . 
which uses fuzzy spheres for compactified dimensions. 

On the other hand, if we wish to obtain a phenomenologically viable 4d model, 
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we may stay in the symmetric phase by keeping the UV cutoff finite and view the 
NC gauge theory as an effective theory of a more fundamental theory. A 0-deformed 
dispersion relation for the photon such as the one displayed in fig. [J] should then 
have implications on observational data from blazars (highly active galactic nuclei), 
which are assumed to emit bursts of photons simultaneously, covering a broad range 
of energy. In particular, a relative delay of these photons depending on the frequency 
could hint at the existence of a NC geometry.l^S 
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